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Quantum theoretical treatment of coherent forward scattering of light in a polarized atomic en- 
semble with an arbitrary angular momentum is developed. We consider coherent forward scattering 
of a weak radiation field interacting with a realistic multi-level atomic transition. Based on the 
concept of an effective Ifamiltonian and on the Heisenberg formalism, we discuss the coupled dy- 
namics of the quantum fluctuations of the polarization Stokes components of propagating light and 
of the collective spin fluctuations of the scattering atoms. We show that in the process of coherent 
forward scattering this dynamics can be described in terms of a polariton-type spin wave created 
in the atomic sample. Our work presents a general example of entangUng process in the system of 
collective quantum states of light and atomic angular momenta, previously considered only for the 
case of spin 1/2 atoms. We use the developed general formalism to test the applicability of spin 1/2 
approximation for modelling the quantum non-demolishing measurement of atoms with a higher 
angular momentum. 
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I. INTRODUCTION 



Various optical phenomena associated with the opti- 
cal pumping process, which have been comprehensively 
studied since the 60's and described in many aspects in a 
famous review [l| by W.Happer, is revived nowadays in 
a new form within the field of quantum information and 
quantum computing. The paramagnetic ground states of 
macroscopic atomic spin subsystems are considered now 
convenient physical objects for mapping and storing the 
quantum information in the quantum states of their col- 
lective angular momenta. The Faraday-type interference 
scheme was proposed for spin squeezing Q and for quan- 
tum communication between atomic ensembles yL The 
proposed ideas were realized in spin squeezing '4' and 
in the entanglement experiments, where a quantum 
measurement on the light forward-scattered from atomic 
ensembles was used. The same kind of off-resonant for- 
ward scattering combined with a quantum feedback was 
used in the recent demonstration of quantum memory 
for light (^]. Applications to various quantum informa- 
tion protocols including a cat state generation and 
quantum cloning of light onto atoms S\ have been pro- 
posed. Efficient generation of entanglement via multi- 
pass interaction have been also proposed Q . Theoretical 
modelling in the above mentioned papers was concerned 
with collective canonical variables for atoms which can 
be conveniently introduced for spin 1/2 atoms. How- 
ever, the actual experiments were conducted using states 
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with higher angular momenta. Hence a theoretical model 
describing off-resonant interaction of light with realistic 
atoms with the angular momentum higher then 1/2 is 
required. 

It is worth noting that the effect of a random Faraday 
rotation due to atomic fluctuations was discussed and ob- 
served first in Refs.fiol llljl more than twenty years ago 
in the context of demonstrating the advantages of light- 
beating method in atomic spectroscopy. The importance 
of atomic polarization was later discussed in Refs. 0, 
where it was shown that the collective spin polarization of 
an atomic ensemble could essentially modify the quantum 
statistics of the outgoing off-resonant probe radiation. 
Quantum statistics of atomic spin variables was first ex- 
perimentally observed via off-resonant forward scattering 
in 14]. 

In the present paper we develop the quantum theory of 
coherent forward scattering on an ensemble of polarized 
atoms with an arbitrary angular momentum. We discuss 
the physical conditions under which the forward scatter- 
ing can be properly described in terms of the effective 
Hamiltonian, generalizing the semiclassical approach P|, 
so that both the field and atomic subsystems are treated 
as quantum objects. Based on the effective Hamiltonian 
we derive the Heisenberg-type equations of motion writ- 
ten for the spatial distribution of the Stokes variables 
of the field and of the collective angular momentum of 
atoms. The solution shows how the entanglement of the 
quantum fluctuations of light and atomic subsystems is 
formed. We show that the entangling process can be 
understood in terms of a polariton-type spin wave in- 
duced in the atomic sample. We support our discussion 
by numerical simulations, assuming the conditions close 
to those of the recent experiments 5] . 
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The paper is organized as follows. In section [H] we 
review the process of coherent forward scattering to 
show how the field Heisenberg operators are transformed 
in interaction with multi-atom ensemble consisting of 
atoms with an arbitrary Zeeman structure and for an off- 
resonant excitation on a dipole allowed optical transition. 
In section Unl we extend the Heisenberg formalism to the 
atomic subsystem and introduce the effective Hamilto- 
nian responsible for the coupled dynamics of an infinite 
number of variables associated with the local polariza- 
tion of the weak probe light and with individual atomic 
spins. The effective Hamiltonian is further transformed 
in section Hvl into another form by introducing the more 
convenient for the polarization sensitive processes irre- 
ducible tensor formalism. In section^we perform meso- 
scopic averaging to arrive at the wave-type equations for 
the collective quantum variables of the light and atomic 
subsystems. These equations are solved and discussed 
in section IVII in the context of the quantum entangling 
process. In particular they are used to test the applicabil- 
ity of the spin one-half approximation for modelling the 
quantum non-demolishing measurement in an ensemble 
of alkali atoms with a higher angular momentum. 



II. COHERENT FORWARD SCATTERING OF 
LIGHT IN THE HEISENBERG PICTURE 

Consider a system of N identical atoms located in a 
finite volume and scattering coherent light of frequency 
uj. The wave of light incident on the atomic ensemble 
is assumed to be weak enough that possible saturation 
effects in its interaction with the atoms are negligible. 
The atomic ensemble, in general, represents an optically 
long (thick in refraction, but thin in absorption) medium 
for multiple off resonance scattering, but, on the average, 
the atoms are separated by distances much larger than 
the wavelength X. Thus each atom is located, on average, 
in the radiation zone of its neighbors. The interaction of 
the atoms with incident and multiply-scattered light is 
assumed to be of the dipole type and a proper description 
of multiple scattering has to be restricted to the rotating 
wave approximation. 

To follow as clearly as possible the analogy between 
classical and quantum descriptions of the coherent for- 
ward scattering it seems convenient to solve the quantum 
problem in the Heisenberg picture. In such an approach 
the basic process is the transformation of the electric field 
operators in a single scattering event. This process is re- 
viewed in appendix^ as an example of the transforma- 
tion of the operator of the free electric field of a moving 
atom. As shown there, the basic result can be written in 
terms of the Heisenberg-type microscopic Maxwell equa- 
tion, where polarization response-operator is given by a 
single point-like scatterer. This allows us to make a sub- 
sequent generalization to the situation of an arbitrary 
number of atoms scattering light coherently in the for- 
ward direction. 



We begin our discussion from the single particle mi- 
croscopic problem, described by equation (|A10ll fsee the 
appendix for details). For a system of N identical 
scatterers randomly located in space it would be not so 
easy to generalize (jA10|) to a multi-particle form allow- 
ing to follow precisely the Heisenberg dynamics. For 
simplicity we restrict the discussion to the case of the 
plain incident wave, as in the experiments with gas cells 
@ . For a transparent medium consisting of N scat- 
terers light is scattered in any direction other than the 
forward direction only incoherently. This means that for 
optically thin and transparent medium the probability to 
observe a photon randomly scattered in any solid angle 
different from the direction of light propagation is just a 
sum of partial probabilities of independent scattering by 
each atom. In the Heisenberg formalism this means that 
the non-forward propagating operator waves created in 
the medium by single or multiple scattering events are 
superposed with random phases to a roughly zero sum. 
The outgoing ffux associated with incoherent scattering 
is proportional to the total number of atoms, but stays 
small because of a negligible value of the off-resonant 
cross section. At the same time for light propagating 
in the forward direction in transparent medium there 
is a very strong coherent enhancement of the scatter- 
ing process. Since the Doppler shift in the forward scat- 
tered modes disappears, see (|A8|) . and the Raman shift, 
caused by Zeeman splitting, does not noticeably change 
the phase of these modes, the transmitted light reveals 
strong coherent superposition of all partial contributions 
associated with single and multiple scattering. For dis- 
cussion of the off-resonant light-atom coherent scattering 
beyond the plain wave approximation see [Tsl . IT^ . 

Let us subdivide the joint dynamics of light and atomic 
subsystems into time increments, during which the num- 
ber of accumulated incoherent scattering events is much 
less than the total number of atoms. This simplifies our 
analysis and allows to consider only the modes coherently 
scattered in the forward direction. In the examples we 
are going to discuss below this type of scattering is addi- 
tionally stimulated by the propagating mode, which has 
the quasi-classical nature. Within the rotating wave ap- 
proximation we introduce a carrier frequency w, coincid- 
ing with the average frequency of propagating light, and 
define an averaged wavenumber k — uj/c associated with 
the full set of the modes propagating in the z-direction. 
Then the spectral bandwidth Aw of the continuum of 
the field modes, contributing to the quantum operator 
expansion, should be chosen to be less than u). The ac- 
tual spectrum of the incident radiation, centered at (D, 
has a narrower bandwidth 5lo, which is much less than 
the detuning from the resonance A = u) — where 
is the resonance frequency of the non-perturbed optical 
transition. We will also assume that |A| is much greater 
than natural linewidth 7„, Doppler shift few, and Raman 
shift oJm'm- We neglect the imaginary shift associated 
with the natural linewidth in the energy denominators. 
The dissipation process associated with the spontaneous 
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decay is neglected in our discussion; the assumptions jus- 
tifying this are discussed in section lVl CI For an analysis 
of the role of spontaneous emission for off-resonant en- 
tangling interaction see 0- 

Under the above mentioned assumpltions the basic ex- 
pression for the polarization of any a-th atom contribut- 
ing to the forward scattering (|A15(I can be rewritten as 
follows 



scattering. That is why the averaging of the operator po- 
larizability (|2.2|) gives only the refraction part of the real 
polarizability tensor of normal Maxwell equations. We 
will further discuss the self-consistency of the dynamical 
approach in section IVI CI 



III. DYNAMICS OF ATOMIC SUBSYSTEM 



P, 



where 



iz,t) 



So 



-Za(t))Ef\z,t) (2.1) 



EE' 



|m')(m|('^)(i) 
(2.2) 

is the polarizability tensor and So is the crossection area 
of the light beam propagating through the medium in 
z-direction. Other notations are specified in appendix 
IXI Expression H2.2|l gives the instantaneous value of the 
polarizability tensor for a-th atom, which depends on its 
exact Heisenberg evolution to the moment t. The Heisen- 
berg operator of the electric field on the right hand side 
of equation (|2.1() is assumed to be dressed in the pertur- 
bations by other atoms located in front of any selected 
a-th atom but non-perturbed by the selected atom itself. 
Here we use Cartesian coordinates in the tensor notation 
and sum over each repeated index. The indices i,j take 
values X OT y. 

Let us consider a thin mesoscopic layer located be- 
tween z and z + Az planes and containing a large num- 
ber of atoms. Such an atomic sub-ensemble will scatter 
incoming field coherently in the forward direction by the 
collective polarization 



N 

p(+)(z,i) = ^p('^'+)(z,t) 

a=l 



(2.3) 



where z coordinate is confined inside the layer (z, z-f Az). 
If the medium is split in a number of thin layers along 
the propagation direction the light beam would be subse- 
quently scattered by this layers in the forward direction. 
Then the "coarse-grained" dynamics of the field opera- 
tors can be described by the following Heisenberg-type 
macroscopic Maxwell equation 



(2.4) 

where the polarization operator P^^\z,t) is subse- 
quently given by expressions (|2.1|l - H2.3|l . 

We conclude this section by the following remark. The 
derived macroscopic Maxwell equation is coupled to the 
corresponding Heisenberg equations governing the dy- 
namics of the atomic subsystem, see below, and it cannot 
be extended up to an arbitrary time t. Its validity is re- 
stricted by ignoring the dissipation process of incoherent 



In this section we discuss how the atomic variables are 
modified via the interaction with the forward propagating 
light. We describe the dynamics of slowly varying ground 
state operators, considering at first a single atom illumi- 
nated with an arbitrary quantized optical field. Then we 
generalize the problem to a macroscopic ensemble and 
introduce atomic collective variables. 



Dynamics of a single atom coupled to 
off-resonant field 



Consider an atom at the origin of the coordinate frame 
at the initial moment of time and drifting in space in such 
a way that during the scattering event it moves much less 
than a wavelength. Let us define an arbitrary dyadic- 
type operator for the ground state of this atom 



T 



\m'){m\ = \jom'){jom\ 



(3.1) 



where by the second equality we specify the atomic state 
more precisely and introduce the atomic quantum num- 
bers: jo and m, m' are the ground state angular momen- 
tum and its projections on the direction of an external 
magnetic field oriented along the Z axis, which in general 
is different from the direction of the light propagation. 
Based on a perturbation theory we can expand the cor- 
responding Heisenberg operator up to the second order 
as follows 



where 



fit) - foit) + f2{t) -f 



To{t) = e^""-''"'|m')(m| 



(3.2) 



(3.3) 



is the operator in the interaction representation and the 
second term 

f2{t) = -\ I dt" [ dt' 
fi Jo Jt" 

X d(t")Eo(i"), [d(OEo(t'),To(t)]](3.4) 

is the second order correction induced by a dipole-type 
optical interaction, see ljA2|l . 

The integral H3.4|l should be evaluated in the rotating 
wave approximation by keeping only the leading terms 
in the limit t ^ These terms can only de- 

pend on normally ordered products of the creation and 
annihilation field operators. In general such products 
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expand over all the spatial modes, but in reality only 
those modes, which will not vanish after the averaging 
over the initial state, are important. These are planar 
modes propagating along the z-axis. Taking into account 
only the propagating modes, the expansion (|3.2(l can be 
rewritten as follows 



fit) » fo(i) + / dt' - Heffit'), Tb(i) 



(3.5) 



where we introduce an effective interaction Hamiltonian 
in the interaction representation 



^(-) 



iz,t)a,,{Lj,t)E^^+\z,t) (3.6) 



Strictly speaking the electric field operators and the op- 
erator of the atomic polarizability tensor should be un- 
derstood here as defined in the interaction representation 
and marked by the index zero. For the electric field these 
operators W^\z, t) are given by the expansion IjAf |l . de- 
fined at z — > 0, where we keep only the forward propa- 
gating modes. The operator of the polarizability tensor 
aij{ijj, t) is given by Ea. H2.2|l . transformed to the interac- 
tion representation, with Q and k = lu/c being the carrier 
frequency and wavenumber of the modes interacting with 
an atom. As in the Ea. (|2.1|l the tensor indices inside the 
definition of the effective Hamiltonian (|3.6|l can run only 
two projections either x or y. 

Note that in our derivation of the effective Hamiltonian 
we used a rather short time increment t consistent with 
the perturbation theory approach; with the assumption 
that the atom does not noticeably change its location 
during the scattering event. Then there would be no dif- 
ference between the interaction and the Heisenberg rep- 
resentations in the evaluation of the integral H3.4|l and in 
introducing the effective Hamiltonian (|3.6|) . Based on a 
general principle of dynamical evolution we can straight- 
forwardly generalize Eq. (|3.()l) up to an arbitrary moment 
in time if we substitute there all the operators in the 
Heisenberg representation. We should also take into ac- 
count a classical drift of the atom in space and consider 
z — z(t) coordinate as its actual location at moment t. 



B. Generalization to a multi-atom ensemble 

Consider now an ensemble consisting of many atoms 
scattering incident light coherently in the forward direc- 
tion. Although the most interesting polariton-type so- 
lutions obtained further in the paper are destroyed by 
atomic motion, we include the motion in our model for 
the sake of generality. However we ignore any possible 
correlations between their spatial motion and the internal 
state evolution. Then each atom is in the environment of 
the field scattered by the atoms located in front of it and 
coupled with such a dressed field via the partial effective 
Hamiltonian (|3.()|) . The full effective Hamiltonian for the 
whole ensemble interacting with the propagating field in 
the Heisenberg representation is given by 

JV 

Tieffit) = - J2 Et\^a{t),t) a[;\0j,t) E^+\za{t),t) 

a=l 

(3.7) 

We preserve here the notation for the full effective Hamil- 
tonian, which was used in the previous equation in the 
case of one atom. 

Then operator H3.f () considered for each atom of the 
ensemble satisfies the following Heisenberg equation 



where a = 1 A^. 

If we take into account the commutation relation be- 
tween Heiesenberg operators of the electric field propa- 
gating in the forward direction, which are given by 



Ei-\z\t), El+\z,t) 



2ttTiuj 



5m5{z'-z) (3.9) 



then the Maxwell equation H2.4|l can be also rewritten in 
terms of the effective Hamiltonian 



J 



dz^' 



1^ 



r 



2u) 



(3.10) 



This operator equation shows that, in turn, the state of 
the field existing at each spatial point of the medium is 
defined by the atomic operators in the Heisenberg repre- 
sentation, according to their spatial distribution in space 
at a given time t. 

Let us make the following remark concerning the defi- 
nition of 5-functions in H3.9|) and in H2.I|I and the validity 
of the dynamical equations in the form H3.I0|) . Our de- 



scription of the quasi-resonant radiation propagating in 
disordered medium in the forward direction is based on 
the rotating wave approximation. According to this ap- 
proximation the spectral bandwidth of the field Aw is as- 
sumed to be at least less than carrier frequency lu. Such a 
truncation of the infinite field continuum makes possible 
to consider the commutation relations of truly Heisenberg 
operators in the form (|3.9|l . Thus the (5-functions here 
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and in H2.1() should be correctly understood as distributed 
in a small mesoscopic area of the order c/Alj. This spa- 
tial scale is obviously longer than A but still much shorter 
than the sample size or than any internal macroscopic 
scale associated with macroscopic susceptibilities of the 
medium. It is also important to think about equations 
(|3.1U|) as on approximation for more general multi-mode 
Heisenberg-Langevin type equation, where the damping 
processes associated with incoherent scattering would be 
taken into consideration. As was mentioned before, the 
vahdity of the output Maxwell equations H2.4|l . H3.10|l is 
actually restricted by the narrow spectral domain 6uj as- 
sociated with the spectrum of the probe light. But in 
the process of transposing the field operators by means 
of the commutation rule (|3.9|) the whole field spectrum 
Auj should be taken into account. 

The coupled equations (|3.8() and H3.10|l considered to- 
gether is the main result of this section. They reveal the 
joint dynamics of the field and atomic subsystems inter- 
acting in the limit of non-saturating off-resonant optical 
excitation. Being an example of Heisenberg equations 
they are valid for arbitrary initial quantum state of the 
field and of the atomic ensemble. The main restriction 
comes from the model of lossless coherent scattering. But 
even with such a simplification these equations are quite 
complicated since they are operator equations for an in- 
finite number of the field and atomic variables. In the 
following sections we identify and discuss special condi- 
tions when equations (|3.8|) and (|3.1U|) could be converted 
to a finite number of Heisenberg-type equations for col- 
lective variables of the field and atomic subsystems. 



IV. REPRESENTATION OF IRREDUCIBLE 
COMPONENTS. STOKES OPERATORS OF THE 
ELECTROMAGNETIC FIELD 

A. Transformation of the effective Hamiltonian to 
the irreducible representation 

The set of dyadic-type operators 1)3. l|l for each atom 
of the ensemble can be replaced with another set of op- 
erators 



(a) 



KQ 



|jom')Oom|(") 




2^ ^jomKQ-^KQ 



(4.1) 



Being the linear combination of original operators 
weighted with Clebsh-Gordan coefficients C " , the pro- 
jectors become irreducible tensor operators, which 
possess the simplest properties with respect to rotational 
transformations, see |l7l |. The linear transformations 
(14. 1|) are consistent with any coordinate frame. Let us 
abandon the choice of the frame with Z-axis along mag- 
netic field, used in the basic definition H3.1|l , and return to 



the original frame with the z-axis along the propagating 
beam, which is more natural for the further discussion of 
the effective Hamiltonian approach. 

For large frequency detuning the Zeeman splitting in 
the energy denominators can be neglected compared with 
the average detuning A — u) — ujjjg , where ujjj^ = loq is 
the transition frequency between the ground and excited 
states characterized by the angular momenta jo and j 
respectively. Then the effective Hamiltonian can be ex- 
pressed as the sum of three terms 

Heffit) = -H%{t) + H« (t) + n%it) (4.2) 

where dependence on time emphasizes that the contribut- 
ing atom and field operators are taken in the Heisenberg 
representation. 

The first term in H4.2|l couples the atomic population of 
the whole Zeeman multiplet and their longitudinal align- 
ment with the full photon flux of propagating light 



HS^(t) 



2Trhu! 



N 

ao(u;)E^oo^W 



5*0 c 



N 



Eo{zait),t) (4.3) 



where 



So(z,i) 



Sqc 
2'Khu) 



E(-)(z,t)E(+)(z,i) 



(4.4) 



is the Stokes operator of the total photon fiux at the 
spatial point z. The isotropic polarizability of an atom 
is given by 



ao(t^) 



Id. 



■303 I 



(4.5) 



where d,-, 



1 ~h{uj-ujjj^) 

is the reduced matrix element of the atomic 



dipole moment for jo — > j transition. The align- 
ment component of the atomic polarizability is given by 
Ea. (|0^ below. 

The second term in Eq. H4.2|l couples the gyrotropic or 
orientation component of the atomic ensemble with the 
Stokes component responsible for circular polarization of 
propagating light 



where 



2'kTiuj 
So c 



N 



,{u^)Y,f[^{t)S,{z,{t),t) (4.6) 



^2{z,t) 



Soc 
2'K'h uj 



a=l 



E^-\z,t)E'^'{z,t) 



Ei-\z,t)Ei+\z,t) (4.7) 



is the Stokes operator associated with circular polariza- 
tion. It is defined in terms of the photon flux operators 
at a spatial point z and shows imbalance between the 
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right-hand (i?) and the left-hand (L) polarizations of the 
field. The oricntational polarizability of an atom is given 
by 



1 



1 1 1 



V2 i Jo ia j j -h{uj-u)jja) 



(4.8) 



As follows from the Maxwell equation (|3.10|l in the clas- 
sical limit, the term Ti-l^fj completely defines the Fara- 
day rotation or other gyrotropy effects existing in a bulk 
medium. 

The third term in Eq. (|4.2|l couples the alignment com- 
ponents of the atomic ensemble with the remaining two 
linear polarized type Stokes components of propagating 
light 



Sqc 

pi' 



a=l 



(4.9) 



where 



2i 



rfiia). 



(4.10) 



The Stokes operator ^3(2, i), showing imbalance between 
the photon fluxes of the modes linearly polarized along 
X and y axes, is given by 



S3(z,i) 



2'Khu] 



-M-Hz,t)E(+\z,t)] (4.11) 



The Stokes operator Si(2;, i), showing imbalance between 
photon fluxes of the modes linearly polarized along ^ and 
r] axis rotated with respect to x and y directions by 7r/4 
angle, is given by 



Si(z,i) 



Sq c 

2'Khu] 



-^(-)(z,i)^(+)(z,i)l (4.12) 



The alignment component of the atomic polarizability is 
defined as follows 



a2{Lu) = (-) 



i-i-jo+i 



\d 



100 1 



1 1 2 



Jo Jo J J -h{ui~^jj,) 



(4.13) 



As follows from the Maxwell equation H3.10|l considered 
in its classical limit the term Ti-pJ^ is responsible for the 
optical birefringence effects with respect to either x and 
y or ^ and 77 directions. 



B. Dynamical equations driven by the effective 
Hamiltonian 



Equation H3.8|l in the irreducible representation trans- 
forms into 



h 

+ 



h 



(4.14) 



where Hq is the unperturbed Hamiltonian responsible for 
the interaction with the external magnetic field. Recall 
here that in general the direction of the magnetic field 
does not coincide with the z-axis. The set of equations 
l|4.14(l for a = 1 TV is equivalent to the set of equa- 
tions l|3.8|l but it is written for more convenient physical 
observables. The irreducible components of a low rank, 
which contribute to the effective Hamiltonian and couple 
to the Stokes observables of the propagating light, al- 
low for a more transparent interpretation than the origi- 
nal projector operators (|3.1|l . The zero rank component 
Tqo ^ (t) is the Heisenberg operator of the total population 
of an a-th atom in its ground state. It can be straightfor- 
wardly verified that the right hand size of H4.14|l is equal 
to zero in this case and the zero rank irreducible operator 
is, in fact, the identity operator. The first rank compo- 
nent T^Q{t), being the Heisenberg operator of atomic ori- 
entation, is equivalent to the vector of the atomic angular 
momentum. This vector undergoes dynamical evolution 
(regular precession and coupling with the field variables) 
caused by the perturbation of the atomic ground state 
by the propagating light. The second rank component 
(t), being the Heisenberg operator of atomic align- 
ment, is equivalent to the ground state quadrupole mo- 
ment of the atom. The quadropole tensor also undergoes 
the dynamical evolution caused by interaction with the 
propagating light. Other higher rank irreducible com- 
ponents are important only as long as their evolution is 
dynamically coupled with the evolution of the lower rank 
components in the complete set of Eas. (|4.14|l . 

The equations H4.14|l have to be solved together with 
the Maxwell equations. Let us introduce the slowly vary- 
ing amplitudes of the Heisenberg field operators 



E\-\z,t) = el{z,t)e'^'-''^' 



(4.15) 



By substituting these expressions into equation (3.10) 
and into its Hermitian conjugated form, we obtain the 
following first order differential equations for slowly vary- 
ing amplitudes 



'd_ ld_ 

dz c dt 

d_ ld_ 

dz c dt 



h{z,t) = — Ueff{t), ii{z,t) 



4iz^t) = ^ [neffit), el(z,i)](4.16) 
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In such Heisenberg-transport equations, as well as in 
the dynamical equations (|4.14|l . we consider time t in 
a coarse- grain temporal scale much longer than |A|~^ 
and the coordinate z on a coarse-grain spatial scale much 
longer than A. Any changes in the atomic subsystem and 
displacement of atoms during the time increments com- 
parable with |A|~^ as well as any changes of the slow 
varying field operators on the scale of a few A are ignored. 
We do not actually know the exact behaviour neither for 
the atom nor for the field Heisenberg operators, but we 
can approximately display their averaged behaviour by 
solution of coupled equations 14.1411 and (|4.1()|) . 

Instead of the Heisenberg-type equations for slowly 
varying amplitudes ei{z,t) and t\{z^t) similar equations 
written for the Stokes operators can be introduced 



d_ 
'd'z 



ld_ 

cdi 



(4.17) 



where « = 0, 1, 2, 3. All the terms appearing on the right 
hand sides of equations H4.14|l and (|4.17|) can be expressed 
via the Stokes operators. Moreover it can be straightfor- 
wardly verified that 5o(z,t) stays unchanged as a func- 
tion of t — z /c, because of the conservation of the number 
of photons in the forward scattering process. In turn this 
means that the first term in Ea. H4.3() . i. e. the isotropic 
component of the effective Hamiltonian, can be omitted 
in (|4.17(l since it commutes with all the Stokes operators. 

As follows from the derived equations the main obsta- 
cle on the way to convert the infinite number of operator 
equations to the system of truncated equations, written 
for the collective atomic variables, coupled with the inte- 
gral or averaged field variables, comes from the spatial de- 
pendence of the interaction process. Such a dependence is 
caused by anisotropic terms in the effective Hamiltonian, 
which lead to spatially varying entanglement of different 
polarization modes of light with atomic spins along the 
propagation path. The spatial profile of the Heisenberg- 
type Stokes operators could be considered as uniform in 
S2 component and as accumulating the collective fluctu- 
ations of atomic spins in complementary Si component 
only in one special case of a pure Faraday effect. In this 
case the anisotropic components of the local susceptibil- 
ities of the medium would be zero on average and would 
exist as fluctuations acting only on the Si Stokes com- 
ponent via random Faraday rotation. This can be true 
if fio\t), f^o\t) and TjikO operators would be zero 
on average along the whole interaction cycle. But even 
if we neglect the repopulation optical pumping mecha- 
nism, coming from incoherent scattering, it would be not 
so easy to show any realistic example of a proper atomic 
transition satisfying this condition. 

Such an example could be jo = 1/2 ^ j — 1/2 or 
jo = 1/2 ^ j = 3/2 atomic transitions. Then there is 
no quadropole moment in the ground state in principle 
and the ensemble consists of the atoms perfectly polar- 
ized (oriented) in the direction orthogonal to the propa- 
gation direction of the probe light. The atom with spin 
one-half in its ground state does exist in reality (^He for 




FIG. 1: Schematic diagram showing the polarization response 
in the probe light, transmitted through the spin oriented 
atomic sample, on random gyrotropy (a) and random bire- 
fringence (b) . Both processes are initiated by transverse fluc- 
tuations of the collective angular momentum. 



example), but it is not a convenient object for the po- 
larization sensitive experiments. In the earlier work on 
entanglement and quantum information protocols with 
atomic continuous variables 0) S S IE S realistic 
atoms were modelled as spin one half systems. One of 
the goals of this paper is to analyze the applicability of 
such approximation. 

Below we derive equations of motion in their general 
form and discuss their possible conversion to the finite 
number of wave-type equations written in terms of spa- 
tially distributed collective variables. 



V. DYNAMICS OF THE SYSTEM IN TERMS 
OF COLLECTIVE VARIABLES 

As the most important practical example, we will fur- 
ther discuss the light propagation through atomic en- 
semble prepared originally in the coherent spin state, see 
figure n We assume that atoms fill the cylindrical vol- 
ume with the cross section Sq coinciding with the cross 
section of the probe light beam. The atomic ensemble 
is located in the homogeneous magnetic field with the 
direction orthogonal to the direction of the propagating 
light. Originally atoms are perfectly oriented along the 
magnetic field in the x direction. Probe light is in a 
coherent state linearly polarized along the x direction. 
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Such a geometry is common, e.g., for experiments aim- 
ing at the quantum state teleportation between field and 
atomic spin subsystems H, IM 1^ - 



where Sijk = ±1 in dependence on order of indices i ^ 
j ^ k, equations (|4.17|) can be rewritten as follows 



A. Field subsystem 

Let us consider first the evolution of the Stokes com- 
ponents. Based on the commutation relation 



^i{z',t),^j{z,t) ^ 2ieijkcd{z - z')^k{z,t) (5.1) 



d_ ld_ 

dz c dt 

d_ ld_ 

dz c dt 

d_ ld_ 

dz c dt 



Ei{z,t) 

^2{z,t) 
S3(^,t) 



N 

aiJ2fiS\t)S{za{t)-z)%{z,t) 

a=l 
N 

a2Y.n'${t)5{zS)-z)t,{z,t) 

a=l 
N 

-a^Y.^[l'{t)5{za{t)~z)^^{z,t) 

a=l 

I 



N 



where we introduced the dimensionless polarizabilities 

(5.3) 



Attuj 

- — a^{uj) 1 = 1,2 
bo c 



These equations are valid for any type of initial condi- 
tions and can be simplified for excitation geometry de- 
scribed in the preamble to this section. 

The coherent forward scattering of light linearly po- 
larized along the x direction does not modify the av- 
erage angular momentum orientation of the atomic en- 



a2j2f^';Ht)Sizait)~z)E2iz,t) 



a=l 
N 



a2j2Tf,it)S{za{t)-z)%{z,t) 



a=l 
N 



a2Y.fl':l{t)5{za{t)-z)%{z,t) (5.2) 



semble. Transmitted light also preserves its mean 100% 
linear polarization. This means that in equations H5.2() 

only Txy\t) and S3(z,t) have non-zero expectation val- 
ues. Other operators exist only as fluctuating quantum 
variables. It is possible to linearize these equations by 
substituting Txy\t) and 23(2, i) by their average values 
and leaving in the right-hand side only the linearized con- 
tribution over quantum fluctuations. Then the evolution 
of the Stokes operators is given by 



J 



d_ ld_ 

dz c dt 

d_ ld_ 

dz c dt 

d_ ld_ 

dz c dt 



N 



E,iz,t) 

Mz,t) 
Mz,t) 



-K2 ^2iz,t) 

K2Si(z, i) 




a^E,J2f^^\t)6(za(t)-z) 

a=l 
N 

a2E,J2T^fit)S{za{t)-z) 



(5.4) 



where S3 = Sq is the average value of the corresponding 
Stokes component, which is approximately unchanged for 
the light beam propagating through the sample. The 
third line in Eas. (|5.4|) just indicates this circumstance. 



The coupling parameter 



N 

«2 = a2Y.f(l>6{za{t)-z) 



(5.5) 
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is responsible for birefringence effects, i.e. for unitary 
transformation of linear polarization (defined in respect 
to ^, rj axes) to circular polarization and vice versa. Here 

fiy^ is the averaged and approximately unchanged value 
of the alignment of a-th atom. For any atom, whose 
angular momentum jo is oriented along the x direction, 
the alignment term is given by 



rp{a) 



[15jo(2jo-l)]^/' 



2[2(jo + l)(2jo + l)(2jo + 3)]i/2 



Txy (5.6) 



The microscopic structure of K2 can be averaged over 
a small mesoscopic interval of Az and substituted into 
Ea. (|5.4() as K2 — a2Txyno, where hq is the linear density 
of atoms (i. e. the number of atoms per unit length). 
From the classical electrodynamics point of view, the 
mesoscopically averaged product K2A is the difference be- 
tween refraction indices for the polarizations along x and 
y directions. 

The most important terms in equations (|5.4|) are the 
last ones on the right hand sides. These terms show how 
any quantum state, originally encoded in the spin fluc- 



tuations of the atomic subsystem, can be mapped onto 
the polarization state of the light subsystem. Thus these 
terms are responsible for the quantum information pro- 
cessing in light - atoms interaction. 



B. Atomic subsystem 

Consider the dynamics of the orientation vector of 
an a-th atom. Applying the commutation rule for ir- 
reducible tensor operators 



rf,(a) rpib) 

^KQt^K'Q' 



K'' 



K+K' + K", { K K' K" 



' Jo Jo Jo 

( \2jo+K" f-,K"Q" f{a) ,r 

see Ref. T?] , to the right hand side of equations (|4.14l) we 
obtain 



^lo^(i) 



+ a2 (-)'■'« VlO 



H„(t), tSwI ± iSj (-)«" — 



2 iLoI h°'^"'"-"^--« 



T «a2 (-) 



± iai 



x/3 



2[jo(jo + l)72jo + l)]^/^ S,(z„(t),t)r(»\(0 



(5., 



As one can see these equations are not closed, since the 
right hand sides are expressed in terms of alignment com- 
ponents. Moreover the higher rank multipoles drive the 
dynamics of the alignment components. Thus to intro- 
duce the system of closed Heisenberg equations it is nec- 
essary to consider the coupled dynamics of all irreducible 
components defined for each atom. 

However the dynamics of the atomic subsystem can 
be approximated by the dynamics of orientation compo- 
nents only as far as we are restricted to the discussion of 
the excitation regime described in the preamble to this 
section. As shown in the appendix IbI equations (|5.8|l can 
be transformed to the set of nonlinear equations, where 
the right-hand sides are expressed in terms of the oper- 



ators of the angular momentum. In a special, but the 
most important for us, regime of small fluctuations such 
transformed equations can be linearized and simplified to 
the equations describing the dynamics of the vector of the 
angular momentum in a closed form. Let us define the 
vector of the collective angular momentum of ensemble 
as 

N 

Then starting from equations (|B7p . written for any a-th 
atom, and making the sum over all atoms of the ensem- 
ble, we arrive at the following equations governing the 
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dynamics of the collective angular momentum 



as follows 



N 



Jz{t) « (r!o + ^2)Jy{t) ~ T,y d2 J2 ^lMt),t) 

1 ^ 

Jy{t) « -{rta + n2)Jz{t) + -T,ai ^ E2izait),t) 

a=l 

« (5.10) 

where parameters ^1,2 and are defined by expressions 
(jB8|l and IjBgp respectively and T^y is given by H5.6|l . The 
last source-type terms on the right hand side of these 
equations are responsible for the mapping of any quan- 
tum state, originally prepared in the Stokes fluctuations 
of the field subsystem, onto long lived atomic spin sub- 
system. 

Equations (|5.4|) and (|5.10|) . considered together, ap- 
proximate the dynamical evolution of coupled collective 
variables of light and atomic subsystems. But as fol- 
lows from the structure of these equations they are still 
not closed because the coupling terms are not expressed 
by collective variables. As we see from (15. 4f) only the 
atoms currently located before the wavefront contribute 
into formation of the fluctuating Stokes variables. In 
turn, such spatially dependent Stokes fluctuations, which 
stay actually unknown, drive the dynamics of atomic col- 
lective angular momentum via source terms in H5.10|l . 
However, as we show below, under certain simplifying 
assumptions, these equations can be further transformed 
into the system of closed equations describing the wave- 
type spatial and temporal distribution of the collective 
Heisenberg operators of the atomic and field subsystems. 



C. Mesoscopic averaging 

If atoms are slowly drifting in space and during the 
interaction with a short probe light pulse each atom pre- 
serves its location inside the area much less than the 
length scale comparable with or with the sample 
size, equations (15. 4() and (|5.10() can be transformed into 
a closed form. Let us note here that if in the experiment 
the duration of a probe pulse is chosen less than a few 
microseconds such a condition is normally fulfilled even 
in case of atoms at room temperature. For the case of 
cold trapped atoms this assumption is consistent with 
the pulse duration up to a second. Then the equations of 
motion can be rewritten for any mesoscopic layer of the 
sample, which gives only a small increment to the Heisen- 
berg operators but contains a large number of atoms. 
The atoms do not leave the layer during the interaction 
time and cooperatively interact with the electromagnetic 
field. If this layer has a length of Az, we can introduce the 
averaged Stokes operators Ei{z,t) along this layer. As a 
next step, instead of total angular momentum, given by 
Eq. (15. 9|) we can define its mesoscopic spatial distribution 



E 

Z<Za<Z-\-Az 



(5.11) 



where the sum over a is extended only over the atoms lo- 
cated inside the layer. Then the total angular momentum 
of the ensemble is expressed as 



3{t) 



J(z, t) dz 



(5.12) 



where dz = Az and L is the length of the sample. 

The above assumptions lead us to the following set of 
closed and coupled equations for the mesoscopically aver- 
aged spatial distributions of the field and atomic variables 



-K2E2{z,t) + 2(3%J,{z,t) 
K.2^i{z,t) - 2eE3 Jy{z,t) 
VLjy{z,t) - 6l.yJ^Si(z,t) 



r d 


1 d' 




az 


^ ~cdt 




' d 


1 d' 




oz 


cdt 





d_ 
dt 



-njziz,t) + 9zX^2{z,t) 

(5.13) 



where we used the same notation for the averaged Stokes 
variables as for their microscopic origins and denoted the 
mesoscopic spatial distributions of angular momentum 
components as J'fj_{z,t), with /i = x,y,z. It is taken 
into account that J^x{z,t) — Jx — const. Equations 
(|5.13(l should be accompanied by corresponding initial 
and boundary conditions, which are given by 



Si(0,t) 

S2(0,t) 

J,(z,0) 



ETit) 



(5.14) 



The solution of these equations presented in the next 
section shows how the swapping of quantum fluctuations 
between light and spin subsystems takes place during the 
interaction process. 

Several new parameters appear in equations l|5.13(l . 
Firstly, by f2 we denote the frequency = f^o + ^^2 of 
the regular precession caused by the external magnetic 
field as well as by the light-induced shift of the Zeeman 
sublevels. Secondly, there are two new parameters in the 
first pair of equations describing the transformation of 
the Stokes variables. The angle (3 



f3 



V3 



2[jo(jo + l)(2jo + l)]i/- 



■ ai 



(5.15) 



is the angle of Faraday-type rotation of the polarization 
plane of the propagating light per one spin flip in the 
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ensemble in z-direction. The parameter e 

^ [15(2jo - 1)]V^ 

' 2[2jo(jo + l)(2jo + l)(2jo + 3)]i/2 



a2 



(5.16) 



is the ehipticity induced in the propagating hght by the 
atomic sample per one spin flip in y-direction, see fig- 
ure ^ Thirdly, there are two angles in the second pair of 
equations describing the dynamics of the spatial distribu- 
tion of the atomic angular momenta. Angle 9y = e is the 
rotation angle of the local collective angular momentum, 
originally oriented along x-axis, around y-axis per one 
photon propagating trough the sample in either ^-type 
or 77-type linear polarization. In turn, angle 0^ = /? is the 
rotation angle of the local angular momentum around z- 
axis per one photon propagating through the sample in 
either right-hand-type or left-hand-type circular polar- 
ization. 



VI. ENTANGLEMENT OF THE QUANTUM 
STATES OF LIGHT AND ATOMS 



For pedagogical purposes, we consider at first a spe- 
cial example of optical excitation in the far-off-resonance 
wing of jo = 1/2 ^ j = 1/2 or jo = 1/2 ^ j = 3/2 
optical transitions and will discuss the general case after 
that. 



A. Example of jo = 1/2 -^j — 1/2, j — 3/2 optical 
transitions 



In the case of jo = 1/2 j — 1/2 or jo = 1/2 
j = 3/2 transition equations H5.13|l are simplified to the 
following form 



2/3S3 J.(z,i) 




" d 


1 d' 




dz 


^ cdt 




' d 


1 d' 




az 


^ cdt 





d_ 
di 
d_ 
di 



Jy{z,t) 



no Jy{z,t) 
-no Jz{z,t) + PXE2{z,t) 

(6.1) 



where it was taken into account that 02 = and therefore 
K2 = 0, 1^2 = 0, e = = 0. Because of the second line 
in the system H6.1|l . one has 'E.2{z,t) — 'E.™{t — z/c). 

The straightforward solution for the distribution of the 
angular momentum components leads to 



Xiz,t) = cosnotjriz) + sin^ioi J^W + 1^^^ I 

Jo 

- sin riot jriz) + cosrioi Jy\^) + I 

Jo 



Jy{z,t) 



t 

dt' sinnoit - t') S^"(t'-z/c) 
dt' cosno(t-t') S^"(t'-z/c) 



(6.2) 



and the Stokes component Si(z,t) is given by 



Si(z,t) = E^\t~z/c) + 2(3Es [ dz' Mz',t-{z~z')/c) = ^^it-z/c) 

Jo 



'2(3% 



dz' cosr!o[t - (2 - z')/c] jriz') + 2/3 S3 / dz' sinr!o[i - [z - z')/c] T^{z') 

Jo 



f-z j*t — {z — z')/c 

2/32 S3 / dz' / dt' sinl7o[i -t' -{z- z')/c] tf{t' - z' /c) 

Jo Jo 



(6.3) 



This solution explains the basic idea of entanglement of 
the collective quantum states of atomic and field subsys- 
tems. After the interaction cycle the quantum state of 
light described by the S™ component is mapped onto the 
atomic angular momentum components because of the 
last terms in Eas. H6.2|l . In turn, the information written 
in either Jl"^ or components of the angular momen- 
tum is mapped onto the Si Stokes component of light, 
because of the second and third terms in Ea. H6.3|l . More- 



r 



over, due to atom- field interaction there is a partial map- 
ping of the S™ component onto Si component because 
of the last term in Ea. (|6.3l) . 

The retardation effects, which are clearly visible in the 
derived solution are mainly important in order to pre- 
serve the proper commutation relations for the whole 
set of the Heisenberg operators involved in the process. 
The initial zero moment of time is a conventional and 
arbitrary chosen parameter to coordinate Schrodinger 
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and Heisenberg pictures. Physically the interaction cy- 
cle starts only when the wavefront of the probe radiation 
crosses z — Q point. We will take this time as the ini- 
tial moment t — 0. In the Heisenberg formalism this 
means that any expectations values of {t) and (t) 
operators or their products are equal to zero for time 
t < 0. Since the goal is to evaluate the expectation val- 
ues of outgoing operators it is acceptable to substitute 
the lower zero limit by t —> z/c in H6.2|l and hy t' —> z' /c 
in the internal integral Ht).3|) . Then the derived solution 
can be further applied at any spatial point z only for time 

I 



lit) = 

Jy{t) = 

and for the output Stokes operator 

sr*(i) = 

+ 2/3S3Cosl)ot + 2/3S3sinl7oi4" 

+ 2(3^%J^ f dt'sinf7o(t-t') (6.5) 
Jo 

which coincides, in principle, with the solution obtained 
earlier in Refs.[a,0|- 

The input-output transformations expressed by (|6.4|) 
and H6.5II have a simple logical structure and therefore 
are attractive for possible applications to the quantum 
information tasks. However, they are only valid for 
jo ~ 1/2 — > j = 1/2, j = 3/2 isolated transitions, and 
it is rather difficult to find a practical example of such 
a transition. However, the same simple input-output re- 
lations are effectively realized in the case of alkali atoms 
excited in the far wing of Di or D2 lines. In this case 
the contributions from the components of the upper hy- 
perfine multiplet to the interaction add up in such a way 
that the alignment effects become unimportant. For a 
transition from the ground hyperfine state Fq to all hy- 
perfine states F we straightforwardly obtain after taking 
the sum over the excited states 

E "2 « E (-)™+i (2F + l)(2Fo + 1) 





1 






j 




{ Fo 


Fo 




{' 


Fo 


Jo J 



Hence the alignment contribution vanishes for the ground 
state electronic angular momentum jo = 1/2 if the fre- 
quency detuning = lu — cupp^ for each hyperfine tran- 
sition Fq ^ F becomes larger than the hyperfine splitting 
of the upper state. The crucial inequality |q;i| ^ \a2\ 



t > z/c, which is in accordance with physical manifesta- 
tion of retardation effects. 

In reality, we are most interested in a correct descrip- 
tion of the behavior of the modes around the carrier fre- 
quency Hj, see definition H4.15|l . For the modes distributed 
in the spectral domain much narrower than c/L, where L 
is the length of the sample, the retardation effects become 
negligible. In this case the above solution is simplified 
and can be written in terms of transverse components of 
the collective angular momentum 



(6.4) 



which allows to neglect the alignment effects is fulfilled 
the better the less important is the dependence on the 
upper hyperfine momentum in the denominator of (|6.6|l . 
In spite of the fact that for heavy alkali atoms the ra- 
tio of the hyperfine to fine splitting in the upper state 
SEhf/SEf is really small (it is less than 10"** for ^^^Cs), 
which makes it possible to neglect the alignment contri- 
bution for large detunings, such an approximation is not 
always easy to fulfill under typical experimental condi- 
tions. This is because for a large frequency detuning one 
has to use pulses with a very large photon number in 
order to make the atom-field interaction strong enough. 
Large photon numbers are often inconvenient due to ei- 
ther limitations of detectors or due to a limited pulse 
duration, or both. Therefore the frequency detuning A 
is frequently chosen to be comparable with 6Ehf. A typ- 
ical detuning chosen in the experiment |^ , is of the order 
of 900 MHz with the hyperfine splitting in &^P^/2 of 
5Ehf ~ 200 MHz. This is one of the motivations for the 
deeper analysis of equations (|5.13|) in their general form. 



B. General solution 



Equations H5.13|l can be solved in the general case, 
which we do first ignoring the retardation effects. This 
is the most important case for practical applications to 
quantum information processing. Without retardation 
the system of coupled equations H5.13|l can be rewritten 



cosl^ot + sinf^ot Jf + /3Jx / dt' sin^a{t ~ t') 3^\t') 

Jo 

- €\^?LQt -f cosl^o^ Jf + / dt' cosl7o(t - t') S™(t') 



r 
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as follows 
az 

a . . _ . 

—E2{z,t) = K2Si(z,i) - 2e^3Jy{z,t) 

d . 

QlJz{z,t) = njy{z,t) - 9yJa,El{z,t) 

-QlJv{z,t) = -nj,{z,t) + 9,J:,E2{z,t) (6.7) 

This is a set of differential equations with constant co- 
efficients and it is accompanied by initial and boundary 
conditions H5.14|l . The solution can be presented in the 
form of an integral response on these conditions. For the 
Stokes operators it can be written as 

%{z,t) = f dt' M^Az,t-t')Ef{t') 

+ f dz' TUz-z',t)Jl''{z') (6.8) 

u=z,y •'0 

with i = 1,2. For the operators of the spatial distribu- 
tions of the angular momentum components the solution 
is 

J,(z,i) ^ Y f dt'g,,(z,t^t')E]p{t') 
3=1,2 

+ J2 f dz'U^,{z-z',t)j:^{z') (6.9) 

with fjL = z,y. The solution can be formally extended 
over the space-time region < i < 00 and < 2: < 00. 

As shown in the appendixlO the kernels of the integral 
transformations in (|6.8|l and H6.9|l can be expressed via 
the corresponding Laplace images 

yZTTl ) J pQ—ioo J sq—ioq 

(6.10) 

where /C(z,t) is any of the matrix functions 
A4{z,t), . . . ,Af{z,t) in equations H6.8|) and H6.9|l and 



JC{p, s) is its Laplace image. The limits of integration po 
and So are chosen as arbitrary real values warranting the 
existence of Laplace images for the spatial and temporal 
transforms respectively. In the appendix [CI the images 
for all the matrices Ai{p, s), . . . ,Af{p, s) are calculated 
explicitly. Evaluation of the integrals H6.10|l can be 
done numerically for any sample when all the external 
parameters are defined. 

Let us briefly explain how the retardation effects can 
be taken into account. To do this the retardation time 
T = t — z/c should be introduced as an independent vari- 
able for each spatial point z at a certain moment of time 
t. Then the basic system H5.13|l . where all the opera- 
tors are considered now as functions of z and r, can be 
converted to the form (|6.7|) but with the derivative over 
retardation time on the left hand side. The initial condi- 
tions should be also modified and written for the initial 
retardation time r = 0, which corresponds now to spa- 
tially dependent moments in the laboratory time t = z/c. 
As before we will assume that the wavefront of the probe 
radiation crosses z — Q point at zero time. In this case 
we need to know the cooperative atom-field dynamics at 
any spatial point z only for time t > z/c. Then we may 
ignore the interaction Hamiltonian for time Q < t < z/ c 
and chose the initial conditions in (|5.14|l in the following 
form 

J,(z,r-0) - cos(r!oz/c) 

+ sin(r!oz/c) j/;{z) = jriz) 

Jy{z,T^O) - -&in{noz/c)j/r{z) 

+ cos(r!oz/c)4;"(z) = 4"'(z) (6.11) 

There are no changes in the boundary conditions in 
(|5.14|) . since at z = the retardation time r coincides 
with t. Thus in the most general case the solution (|6.8|) - 
(|6.10(l should be rewritten in terms of the retardation 
time and should depend on jf™'(z) and J™'{z) as on 
initial quantum fluctuations of the atomic angular mo- 
menta. 

Finally the retardation effects give us the following cor- 
rection to the spatial and temporal dynamics of the field 
variables 



E,{z,t) ^ Y / dt'M^j{z,t-z/c~t')EJ'{t') + Y / dz'J^,,{z~z',t-z/c)X"'{z') (6.12) 

J = l,2 i^=z,y Jo 





and of the atomic variables 

f^t — z/c 



pt — Z I C (' z 

J^{z,t) ^ dt'g^,{z,t-z/c-t')Ef{t') + Y dz'M^.i 

3 = 1,2 ■'0 u=z,y •'0 



z-z' ,t-z/c)TJ"[z') (6.13) 



These expressions can be applied for calculation of ex- pectation values of any products of operators at a spatial 
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point z for the time t> zjc. 

An important feature of solutions (|6.12l) and Ht).13|l is 
their wave nature. As follows from the Laplace trans- 
form carried out in the Appendix [CI the eigenmodes of 
spatial and temporal distributions of atoms-field collec- 
tive fluctuations are associated with the poles of determi- 
nant A(p, s) 0, where A(p, s) is defined by l|C3|l . For 
atoms with spin 1/2 this condition leads to p ^ and 
s zbirio- In the absence of the magnetic field JIq ^ 
there are only two collective modes associated with the 
total number of atoms and with the total number of scat- 
tering photons which undergo polarization sensitive in- 
teraction in the process of coherent forward scattering. 
This makes it possible to further simplify the effective 
Hamiltonian and to discuss it in terms of integrated col- 
lective variables of the canonical type |^ ■ 

However in a general case, due to presence of align- 
ment associated effects, the condition A(p, s) ^ gives 
coupled roots s = s(p), i.e the wave-type modes which 
can survive in the sample and be excited by either time 
dependent fluctuations of the field Stokes components or 
by spatial fluctuations of atomic spins. Such quantum 
superposition of the field and material waves can be un- 
derstood as a polariton-type spin wave generated in the 
medium. The multimode structure of such waves does 
not allow entanglement only between integrated collec- 
tive variables. Instead, in general the entanglement is 
distributed in all polariton modes. From the quantum 
information point of view this means that in the case of 
atoms with the ground state angular momentum equal or 
larger than 1 a possible entanglement resource hidden in 
the polariton modes allows for multimode entanglement 
rather than a single mode entanglement available for spin 
one half systems. The quantum information then will be 
mapped on spatially dependent correlations of an atomic 
standing spin wave. 

The following remark concerning the terminology is 
due here. Let us point out the principle difference be- 
tween the spin- wave dynamics described by Eqs.(5.13) 
and the collective spin behavior in the spin-organized 
cold fermion gas, which was discussed years ago, see pof. 
In our case the cooperative dynamics of atomic spins is 
driven by the radiation-type interaction whereas in the 
cold fermionic systems it is driven by the static interpar- 
ticle interaction via the longitudinal electric field. 



C. Spin squeezing in an ensemble of cesium atoms. 

In this section we concentrate on calculations relevant 
for the existing experimental example carried out with 
ensemble(s) of -'^^'^Cs atoms. We first calculate the de- 
gree of spin squeezing ignoring the alignment, i.e., follow- 
ing the model adopted in Q and describing the atomic 
ground state by the spin (orientation) only. We then 
present the numerical results including the alignment and 
show how this affects spin squeezing of atoms. 

In the experiment l5l| the entanglement (two mode 



squeezing) was generated between two spatially sepa- 
rated via a Faraday- type detection of light. For peda- 
gogical reasons we will make numerical simulations for 
the case of single mode spin squeezing, which makes no 
difference for the present discussion. We will completely 
ignore the retardation effects and consider the spin dy- 
namics without external magnetic field. 

In the case of cesium atoms the alignment contribution 
can be suppressed if the frequency detuning of the probe 
light from the atomic resonance is much larger than the 
hyperfine splitting in the upper state, as described in 
the end of section IVI Al If the alignment is ignored the 
basic equations stem from the expressions (|6.4f) and (|6.5|l . 
which in the absence of magnetic field read 

JM = JT 

Jo 

sr*(t) = s™(t) + 2PE3 JT 

Sr*(t) = S™(t) (6.14) 

The input-output transformations H6.14|l show the en- 
tangling mechanism of the atoms-field variables in the 
process of coherent forward scattering. If the number 
of atoms and photons is large enough such as pjx S> 1 
and P'E.st ^ 1 the output quantum fluctuations become 
strongly entangled. Indeed, in this case the role of J™ 

and {t) terms on the right hand side of (|6.14() becomes 
negligible if these fluctuations have been originally Pois- 
sonian. After the interaction there would be strong corre- 
lations between fluctuations of J* 'E1^*{t')dt' and Jz{t) as 
well as between Jy{t) and J* 'E.'^^{t')dt' . If the observable 

/p E°^*{t')dt' is measured by a balanced Faraday detec- 
tor there will be no longer standard quantum uncertainty 
in Jz{t + 0), whose original (at t ^ 0) Poissonian fluc- 
tuations will be suppressed by a factor /JSai ^ 1. The 
collective spin state will become squeezed and will accu- 
mulate its quantum uncertainty in Jy{t + 0) fluctuation. 

It is important to recognize that even in such an ideal 
scheme there is a restriction on the number of participat- 
ing atoms and photons because of accumulation of events 
of the incoherent scattering. It is clear that in each event 
of the incoherent scattering one photon and one atom of 
the ensemble cancel out from entangling process. These 
losses can be neglected if the number of scattering events 
is much less than the total numbers of atoms Na and 
incoming photons Nph- This can be written as two in- 
equalities 

^NphNa « Na 
Do 

^NphNa « Nph (6.15) 

DO 

where is the cross section for off-resonant incoher- 
ent scattering in full solid angle, and Sq is the area of 
the light beam, which coincides with the area of atomic 
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cloud. Both inequalities lead to similar restrictions on 
the number of participating atoms and photons 



Na,Nph < 



(6.16) 



In case of photons this inequality can be also under- 
stood as a restriction on the whole interaction time T 
(the probe pulse duration), since Nph = SpT = SsT. 
Since both type of the losses are undesirable one can as- 
sume that Na = Nph to exclude any preference for atoms 
or photons. The number of scattered photons can be 
expressed as 



N, 



ph 



V- 



(6.17) 



where by small parameter rj we denoted the relative num- 
ber of atoms and photons lost as a result of incoherent 
scattering, see the first line of H6.15|l . A detailed analysis 
of the role of spontaneous emission is presented in [9j . 

The scheme of energy levels of ^^'^Cs is shown in fig- 
ure [3 Skipping technical details of calculation of off- 
resonant cross section cta (on an atom in the Zeeman 
sublevel Fq = 4, Mq = 4 of its ground state) and of es- 
timation of the elementary Faraday angle /3, we present 
the final result. The square variance of any output Stokes 
component where i = 1, 2, is given by 
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FIG. 2: Energy level diagram of D2-lme of "^Cs: A5, A4, 
and A3 are the frequency detunings of the probe light from 
hyperfine transitions participating in the process. 



function /(...) denotes the product 



%T [1 + C.(J)] (6.18) 



where the Mandel parameter ^i(J), considered as func- 
tion of the total collective angular momentum J « J^^, 
shows the relative deviation from the shot-noise level. 
This deviation is due to the variance of the atomic state 
mapped onto the output light. For SJ"* Stokes com- 
ponent the original value of the Mandel parameter ^1 
transforms as 



(6.19) 



where the second line defines the dimensionless parame- 
ter K responsible for coupling of the field and atomic sub- 
systems described in terms of their canonical variable, see 
Q • For S™* component the Mandel parameter preserves 
its magnitude, so that ^2{J) = In (|6.19|l and through- 
out we approximate the input fluctuations as Markovian 
delta-correlated process in the lower frequency domain 



(6.20) 



and assume the Poissonian (coherent state) square vari- 
ance of the angular momentum (J^) — Jx/^ — J/2 for 
original coherent atomic spin state. Second term in (|6.19|) 
denotes the contribution of the Faraday effect itself and 
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(6.21) 



where 7 is the rate of the spontaneous decay of the upper 
state and Aj = lo — L0j4 with j = 3,4,5 are the frequency 
detunings of the probe light from each exciting hyperfine 
transition, see figure [3 

Note that, as a consequence of neglecting the align- 
ment effects, the total angular momentum (or a number 
of atoms) contributes in the output Mandel parameter 
only in combination (3J. This parameter is the optical 
activity of the sample, i.e. the angle of Faraday rotation 
of the planar wave for the ensemble perfectly oriented 
along the probe beam, see figure ^ This linearity (see 
dashed-dotted lines in figures has been used as the 
benchmark for the determination of the projection noise 
level of atoms in all the work related to quantum state 
generation with atomic ensembles, see for example 

As follows from the above discussion of the role of 
spontaneous emission, a reasonable choice of parame- 
ters for spin squeezing corresponds to PJ Forif{. . .) = 
4?7/(. . .). Under typical conditions for the frequency de- 
tuning A5 ~ 1000 MHz and 77 '--^ 0.1 the optical activity 
can be PJ 1, which can provide spin squeezing far 
below the standard quantum limit in (|6.18() . 

In order to calculate the corrections to the variances 
of the output Stokes components for a realistic alkali 
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atom at a finite detuning the alignment associated ef- 
fects have to be included in calculations. This can be 
done only numerically by means of the transformations 
(|6.8|) . In figures I3I5I we plot the results of these calcu- 
lations. The variances of the Stokes components S°"* 
and S™* are shown as a function of the optical activ- 
ity of the sample /3J, calculated for different frequency 
offsets A5 = 700, 1000, 1200 MHz and for rj = 0.1. In 
the figures the solid curves represent the complete result 
within the described model, whereas dashed curves are 
computed only including the partial contribution com- 
ing from the first term in Ea. H6.8|l . These curves show 
the excessive atomic spin fluctuations, mapped onto the 
Stokes collective variable of the transmitted light, beyond 
the level of the transformed original light fluctuations 
and beyond the level of the spin (orientation) projection 
noise. The incoming light is taken to be in a coherent 
state with Poissonian statistics and with ^1 = ^2 = 0. 
The enhancement of the light noise is caused by higher 
orders of the interaction process. Considered at certain 
value of PJ the difference between the complete variance 
and the background level shows the efficiency of the Fara- 
day detector as a spin squeezer. The scheme works the 
better the higher this difference is. However, because 
of the alignment the squeezing has to be considered not 
for a collective angular momentum but for a more com- 
plex type of a spatial spin mode, which is defined by the 
input-output transformations H6.8|l and H6.9I) . Squeezing 
and entanglement for such spatial modes will be consid- 
ered elsewhere. 

It is instructive to compare the result of numerical cal- 
culations with the approximation ignoring the alignment 
effects. The dashed-dotted linear dependencies in figures 
I3I5I reproduces the ideal result described by Eq. (|6.19|l . 
As we see, for small values of (3 J the Faraday approxima- 
tion fits the exact solution with rather good accuracy and 
the deviation of ^2 from zero is quite small. This means 
that for the sample characterized by rather small optical 
activity the spin one-half approximation for the multi- 
level cesium atom is self-consistent and potentially good 
for describing the real experimental situation. Hence for 
such a sample the squeezed spin standing wave can be ap- 
proximated by the collective angular momentum. Let us 
point out here that in the existing experiment .5] the op- 
tical activity was (3 J < 0.2 and the alignment correction 
was not so important. However, as clearly seen for large 
f3J the difference between the input-output transforma- 
tions in their general form and approximation (|6.14() be- 
comes quite important. 



VII. CONCLUSION 

We have considered the quantum theory of coherent 
forward scattering of light by an ensemble of multi-level 
atoms polarized in their angular momenta. As a result of 
such a scattering process the quantum states of the field 
and atomic subsystems are transformed into an entangled 
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FIG. 3: The Mandel parameters for the output square vari- 
ances of the Stokes components 2°"' (upper plot) and H2"* 
(lower plot) as a function of optical activity of the sample 
for the frequency offset A5 — 700 MHz. The dotted curve 
indicates the original shot-noise level, the dashed curve is 
the atomic signal coming from the first term in the 
dashed-dotted linear dependence is the Faraday approxima- 
tion lb. 191 . the solid curve is the complete output variance. 

state. In our discussion of the process we followed the ef- 
fective Hamiltonian approach, which in the semiclassical 
form is normally applied for studying optical pumping 
processes after adiabatic elimination of the excited state. 
Compared to earlier studies of entanglement of light and 
atomic ensembles, we have derived the effective Hamil- 
tonian in a more general form for the atoms with an 
arbitrary total magnetic momentum. Towards this end, 
we have found the physical conditions under which the 
analysis can be simplified by introducing a finite number 
of collective variables for light and atoms. 

We showed that under certain conditions the cooper- 
ative atoms-field dynamics can be properly described by 
the wave-type coupled equations for the space-time evo- 
lution of the collective Heisenberg operators of the field 
and atomic subsystems. In these equations an infinite set 
of atom-field operators is truncated via introduction of 
slowly varying collective modes. The coupled equations 
for the time evolution of spatially dependent operators 
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FIG. 4: Same as in figure |31 for the frequency offset A5 — 
1000 MHz. 
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FIG. 5: Same as in figure |31 for the frequency offset A5 = 
1200 MHz. 



of the Stokes components of light and of the macroscopic 
fluctuations of the coUective atomic angular momentum 
are written in a closed form. In the general case the 
coupled dynamics of the atom and field operators man- 
ifests itself in a spin polariton wave created in the sam- 
ple. The fluctuating components of atoms and field be- 
come strongly entangled in the polariton wave in space 
and time. Such spin polariton waves initiated by radia- 
tive forces are different from the collective spi n dy namics 
existing in the spin-polarized quantum gas [23|- They 
are also different from the polariton modes discussed in 
PH where polaritons of the combined atom-light state 
are introduced. In our case the quantum entanglement 
arises from the interaction between the internal collec- 
tive polarization degrees of freedom of light and atomic 
subsystems. In particular our analysis yields the input- 
output transformations for the Heisenberg operators of 
the collective variables of light and atomic spins after 
the whole interaction cycle. Our results suggest that a 
successful approach of using off-resonant light-atomic en- 
semble interaction for quantum information processing 
can become even more fruitful with the use of a multi- 
mode type of entanglement provided by spin polariton 
atomic variables. 



Numerical simulations demonstrate the importance of 
the developed formalism for application to a realistic ex- 
perimental situation. We considered a well known ex- 
ample, when the Faraday rotation is used as a non- 
demolishing measurement of transverse fluctuations of 
the collective atomic spin and can be utilized as a phys- 
ical mechanism for the spin squeezing in an ensemble 
of spin polarized atoms. We tested the validity of the 
spin one-half approximation for a realistic Cesium atom, 
which is normally used to describe the interaction with a 
far off-resonant probe light around Di or D2 transitions 
of alkali atoms. As shown, for small values of optical ac- 
tivity this approximation is self-consistent and deviates 
negligibly from the calculations based on a general solu- 
tion. However for the samples with high optical activity, 
where quantum correlations become strong, there is an 
important quantitative as well as qualitative difference 
between the general solution and the model of spin one- 
half atoms. 
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APPENDIX A: TRANSFORMATION OF THE 
ELECTRIC FIELD IN A SINGLE SCATTERING 

The unperturbed electric field operator in the origin 
of the coordinate frame coupled with a scattering atom, 
which drifts with velocity v, is given by 



V 



1/2 



3-i(wfc-kv)t 



(Al) 



where , aj^^^ are respectively the operators of creation 
and annihilation of the photon with wavevector k and in 
the polarization state ek^. V is the quantization volume. 
The second line in (jAlj) defines negative and positive fre- 
quency components of the electric field. 

The dipole-type interaction operator of the atom with 
the electric field is given by 



(A2) 



where the operator for the atomic dipole moment d(t) 
and for the electric field are defined in the interaction 
representation. Based on a perturbation theoretic ap- 
proach, the exact solution for the Heisenberg operator 
can be written as the following expansion 



E(t) =EoW + Ei(t) + E2(t) 



(A3) 



which can be also written for positive E(+)(t) and nega- 
tive E*^~^(i) frequency components. We assume that the 
wave function of the joint atom-field system describes the 
combined state where an atom occupies the ground state 
and the electric field is in a weak quasi-coherent state 
non-saturating the atomic transition. Then the correc- 
tion of the first order in ljA3|l will disappear after av- 
eraging over the wave function. Thus the second order 
correction gives us the main contribution, since it is re- 
sponsible for the scattering process. The second order 



dt" 



dt' 



X d(i")Eo(i"), [d(t')Eo(i'), Eo(t)JJ(A4) 
As follows from this expression, in a complete dynami- 
cal description of the process there is a memory of initial 
conditions in the formal expansion of perturbation the- 
ory. However, for non-saturating fields this solution can 
be spread out over the time t 3> 7~^, where 7 is the nat- 
ural radiative relaxation rate of the upper state. But in 
this case it is necessary to take into account all the ra- 
diative correction for the retarded and advanced Green 
functions of the decaying upper atomic state. This can 
be done by introducing the natural decay law into the 
time behavior of these functions. Then the integral (|A4|I 
loses its dependence on the lower limit as on initial time 
coordinate. Let us also point out that there is only con- 
ventional choice of initial time <o ~^ to coordinate the 
Heisenberg and Schrodinger representations. The real 
physical conditions can be arranged as the wave front of 
probe radiation (expressed by expectation values of any 
products of freely Heisenberg operators of the field) could 
arrive to the interaction area at any time after to = 0. 
Thus we can always think that observation time t satisfies 
inequality t 3> 7"^. At the same time it is important to 
recognize that by including the radiation decay into the 
Green functions of the upper atomic states we average 
their time evolution and ignore any random fast varia- 
tion of the field and atomic operators associated with 
high frequencies of the field continuum. 

The integral ljA4p gives us the solution for the electric 
field in the origin of the frame coupled with a moving 
atom. But in the zero order of relativistic effects, when 
only retardation effects in the radiation zone have to be 
taken into consideration, this solution coincides with the 
electric field in the laboratory frame at the point of atom 
location as well as in the small vicinity of this point. Then 
we can obtain the solution for any point in the labora- 
tory frame by using the propagation law in free space. 
By this procedure one obtains the following expansion 
for the positive frequency component of the electric field 
operator in the radiation zone (r ^ A) of the scattering 
atom 



E(+)(r,i) =E[,+'(r,t) + E^+'(r,t) + . 



(A5) 



where 
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E(+)(r,<) 



[ V 



(A6) 



and 



E(+)(r,i) 



EEE 
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m')(m| flk^ 



i{d±)nm (dek^)m'„ i(d^),„'„ (dek^) 



i{ui' + ujnm - k'v) i(cjfe - ti;„„i - kv) - 7„/2_ 



(A7) 



and the negative frequency component is given by the 
Hermitian conjugation: E'^~^(r,t) = E(+^^'(r,t). Here 
the origin of the frame is chosen in the location of the 
atom, which is assumed to be unchanged during the hght 
propagation time r/c. The scattered hght frequency uj' is 
defined here via the input frequency Uk , and the Raman 
shift uJm'm for atomic transition |m) |to') There is an 
additional Doppler shift caused by atomic motion, given 
by 



LU' = UJk' 



(k' - k)v 



(A8) 



The wave vector k' = uj'r/cr, but in the right hand side 
of l)A8|l the Doppler correction is neglected and it is as- 
sumed that k' « Lor/cr. The transition dipole moments 
in l|A7p are defined in the Schrodinger representation and 
its transverse component is given by 



(A9) 



The sum over n is expanded over all possible excited 
transitions characterized by natural linewidths 7„, but as 
a practical matter, the sum can be restricted to the most 
significant resonance transitions and the frequency cuk = 
u can be associated with the frequency of the incident 
mode. 

In spite both the terms in ljA6p are well known in the 
scattering theory, see J^], we claim that only the lat- 
ter term, constituted with rotating wave approximation, 
should be left in aiming the following generalization to 
multiple scattering process. Indeed, the fist term con- 
tributed in (|A6|) reveals the quantum mechanical phe- 
nomenon than the scattered photon is created in ad- 
vance than incoming photon was annihilated. This pro- 
cess has very small, and in fact negligible, amplitude and 
should be ignored if one stays with semiclassical-type un- 
derstanding of multiple scattering process as sequence of 
successive scattering events. Moreover, as shown in [23|. 
the electric field transformation in the presence of one 
scatterer in form l|A5p - l|A7|l preserves commutation rela- 
tion, i. e. being unitary transformation only in rotating 
wave approximation. 

In the rotating wave approximation the perturbation 
theory solution IjASp - ljATp satisfies the following equation 



AE(+)(r,t) - l|!,E(+)(r,i) 



,PV+'(r,t) 



where 



(AlO) 



pW(r,t) 



(27r)3 



and in turn 

pi-)(k',t) = EEE 

-icj t I™/ 



V 
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X e 



to')(to| iflk^i 
{d±)m'n (dek^) r. 
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(A12) 



is the operator of transverse component of atomic po- 
larization responding on external field. Frequency u' is 
given by Eq. ljA8|) , but wavevector k' is independent vari- 
able here. 

In above discussion the interaction time t was assumed 
to be short enough for validity of the perturbation theory 
approach. Therefore in the derived equations the dyadic- 
type operators |m')(m| of the low atomic state can be 
selected as referred to interaction representation 

|m')(m|o(i) = e''^™'"* |m')(m| (A13) 

But this operators only slightly differs from their Heisen- 
berg analogs during the short interaction time. There- 
fore the evolution of the electric field operator can be 
extended up to arbitrary time i, if the final equations 
of this section are modified as follows. All the dyadic- 
type operators should be changed by the corresponding 
Heisenberg operators with keeping the complete dynamic 
evolution up to moment t 



|m')(m|o(i) ^ |m')(m|(t) 



(A14) 



The drift of the atom in space cannot be further ignored 
and instead of origin of the laboratory coordinate frame 
we should assume its actual location associated with its 
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classical motion Va = J^a{t). Equation (jA10|) stays valid 
and unchanged but atomic polarization HA12I) modifies 
to 



771, m' 71 k/j- 

(d±)m'n (dek^) 

nm 



2'KhLOl 
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(A15) 



and Va = Ta (t) ■ In substituting HA15p into Eq. (|A10p the 
time derivation should be done only for the fast oscillat- 
ing components of the Heisenberg operators and for the 
exponential factors. 



and the alignment tensor is given by 



f (a) ^ Vl5 

[2jo(jo + l)(2jo - l)(2jo + l)(2jo + 3)] 1/2 

^ "-^Iglg' [Jg Jq' + 

qq' 



where jg"'' are the cyclic components of the operator vec- 
tor of angular momentum, which are defined by their 
Cartesian components as follows 



(B2) 



APPENDIX B: LINEARIZED DYNAMICS OF 
ATOMIC ANGULAR MOMENTA 

The operators of irreducible components can be ex- 
pressed in terms of the operators of atomic angular mo- 
menta. For any a-th atom of ensemble the orientation 
vector is given by 



p{a) 
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(B3) 



By substituting subsequently (|B1|) - (|B3|1 into equations 
(|5.8|) the latter can be straightforwardly transformed to 
the set of nonlinear equations containing only the oper- 
ators of atomic angular momenta 
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(B4) 



In their general form these equations are quite compli- cated and not closed because of their nonlinear structure, 
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but they can be simpHfied and Hnearized in the follow- 
ing assumptions. The dynamics of operators jl""* (t) is 
driven only by those terms, which, being averaged, have 
a quadratic scale over fluctuations of the field and atomic 
variables. Recall that Sq = 53. Physically this means 
that there is no coherent process demolishing the origi- 
nal spin orientation of single atom along x direction. So 
far we neglected any possibilities of incoherent scatter- 
ing our analysis has to be restricted by assumption that 

in average ji°'\t) = jo- Moreover, since this observable 
has a maximal possible expectation value and it has no 
deviation from jo in the lower orders of weak external 
perturbations, it is allowed to approximate the Heisen- 
berg operator jl"^ (t) when it appears in linear or in the 
squared nonlinear form by its non-perturbed projector 
onto atomic wavefunction 

(i) jo I Jo , Jo) (jo , Jo I (t) ~ const* (B5) 

But while substituting it in the operators' products one 
has to follow the rule 

jtHOj^'W +3tH0it)W - (2jo-i)itHt) 

jtHOit^W +3tH0ji"^(i) - (2jo-l)jtHt) (B6) 

This is crucial point of linearizing procedure, it is 
expected that the left and right-hand side operators 
are coordinated if the spin subsystem is slightly dis- 
turbed and the time dynamics of off-diagonal projectors 
|jo,Jo)(jo, jo-l|(i) and |jo, jo - l)(jo, jo|(i) is only taken 
into consideration. Then orientation dynamics, described 
by (jB4p , can be approximated by the following set of lin- 
earized equations 

t\t) « (^^o + n2)jl;\t) - Ti;)a2Si(z„(i),t) 



« 

where 



[15(2jo-l)]i/2 



[2jo(jo + l)(2jo + l)(2jo + 3)]i/2 



(B7) 



as So (B8) 



is the light induced shift between m = jo and m — jo — I 
sublevels. Average value of alignment component Txy\ 
is defined in the body of the paper by Ea. (|5.6|l . Simi- 
larly we introduced here the orientation component Ti"-* , 
associated with average orientation of a-th atom, as 

[3jo]^/^ 



rp{a) 



[(jo + l)(2jo + l)]i/2 



(B9) 



Here f^"^ is the average orientation component of 10- type 
defined in the frame with Z-axis along magnetic field, 
which coincides with x-axis in our case. This component 



is the same for all the atoms of ensemble. Making the 
sum over all partial equations HB7|I we come to equations 
(15.10(1 written for the collective vector of atomic angular 
momentum. 

APPENDIX C: LAPLACE TRANSFORM OF THE 
ATOMS-FIELD DYNAMICAL EQUATIONS 

Let us define the Laplace images of the space-time de- 
pendent Stokes components of the probe light and of the 
collective angular momentum of atoms 

POC POO 

E,{p,s) ^ dzdte-P'-''E,{z,t) i = l,2 

Jo Jo 

/>oc poq 

JMs) = dzdte-P'~''Mz,t) fi = z,y (CI) 

Jo Jo 

where parameters p, s > 0. Then the original system of 
differential equations ((6.7(1 with initial and boundary con- 
ditions ((5.14(1 can be transformed to the following system 
of linear algebraic equations for the Laplace images 

pSi(p,s) + K2E2{p,s) - 2p%X{p.s) = Sl"(s) 
-K2Si(p,s) +pS2(p,s) + 2e^:iJy{p,s) = S^"(s) 

J,Si(p,s) + sj,{p,s) - njy{p,s) = j-^b) 
-ezJ.^2{p,s) + nX{p,s) + s jy{p,s) = j;^{p) 

(C2) 

The determinant of this system is given by 

A = A(p,s) - {s^ + n?){p^ + Kl) 

+ 2{ee, + p9y)J^ S3 sp - 2{eey + pez)Jx S3 n K2 
+ A(3ee,ey{J,^^f (C3) 

and its solution can be easy found by means of the 
Kramer's rules. 

The solution is expressed as linear transformations of 
the Laplace images of the boundary Stokes operators for 
the field and of the initial angular momentum operators 
for the atoms 

%{p,s) = M,M^)^T(') 

i=i,2 

+ J2 ^Up,s)j:^ip) 

Mp^s) = ^ g^,ip,s)E^is) 
i=i,2 

-t- J2 ^,u{p,s)Jl-{p) (C4) 

These transformations perform the Laplace images of the 
integral transforms ((6.8(1 and ((6.9() introduced in the body 
of the paper. 

The self-transformation matrix Jvi (p, s) of the Stokes 
operators is given by 
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M{p, s) 



1 f Pis'' + n') + 2e6, X S3 s, -K2{s' + n') + 2(39, X S3 n 



A(p, s) V «2(s' + f^') - 2e0y X S3 n, p{s' + n') + 2fiey J, S3 s 

I 



(C5) 



The self-transformation matrix A/'(p, s) of the angular momentum operators is given by 
I 



Af(p,s) 



1 



A(p, s) 



2ed,XE3p, n{p' 
2f39,XE3K2, s{p' 



2e0y Jx S3 K2 
2(3eyJxE3p 



r 



(C6) 



The cross-transformation matrix J-'{p,s), which is re- 
sponsible for mapping the initial angular momentum fluc- 



tuations onto the outgoing Stokes components of the 
transmitted light, is given by 



J 



Hp, s) 



2E3{pnp + eK2s) 

A{p, s) V 2E3{enp + I3k2 s), -2S3(eps - ^^2^) - AfieOyX Si 



1 / 2S3(/3p s_- eK2n) + AM.Jx Si 



r 



(C7) 



The cross-transformation matrix Q{p^ s), which is rcspon- spatially distributed angular momentum fluctuations, is 
sible for mapping the input Stokes components onto the given by 

I 



Q(p, s) = 



A{p, s) 



-Jcci0yPs~e,H2n) - 2e9ye,J^E3, 



Jx{0yK2 S + 9z^p) 



Jx{o,K2 s + eyQp), JxiOzPs- 9yK2^) + 2peye,j^E3 



r 



(C8) 



Thus the Laplace images of the field and atomic variables 
become fully defined. 

To return into original space-time dependent represen- 
tation it is necessary to evaluate the integrals (|6.1Q|I . As 
we see from (|C3(l - (|C8|l the Laplace images of all the ma- 



trix elements have polynomial structure and the return 
transform could be found for any definite set of the ex- 
ternal parameters, such as K2, O, . . . etc.. 
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